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Emission Mechanisms in High Energy Pulsars: From 
Gamma Rays to Infrared 
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Abstract. We present a model for the gamma-ray emis- 
sion of pulsars in terms of curvature radiation by highly 
relativistic particles. It is shown that the injection of 
power-law primary particles from an outer gap and sub- 
sequent cooling by curvature radiation losses reproduces 
the high energy spectrum and the luminosity of pulsars. 
As a result a spectral break is to be expected at a pho- 
ton energy of ~ 1 GeV. This value does not depend on 
the surface magnetic field strength or on the period of the 
neutron star, but only involves a geometric factor. The 
predicted change of the spectral index by 1 explains in 
particular the spectral shape of PSR B1706-44. We find 
that according to this model the luminosity of pulsars in 
the high energy band varies according to L 1 oc Bq/P 5 ^ 3 , 
which is in good agreement with the observations. 

A model for the infrared, optical and soft X-ray emis- 
sion of pulsar is also presented. It is based on anisotropic 
synchrotron emission by relativistic particles in an outer 
gap scenario with a single energy distribution AT (7) oc 7~ 2 . 
It is shown that this synchrotron model is able to repro- 
duce the spectral shape between the infrared and soft X- 
rays and also the corresponding luminosities for the Crab 
pulsar. In particular, the long standing problem of the 
rapid spectral decline towards infrared frequencies is un- 
derstandable as emission at very small pitch angles from 
low energy particles with 7 <; 10. It is also shown that the 
scaling of the synchrotron model explains the observed 
correlation between the X-ray luminosity and the spin- 
down luminosity of the neutron star L x ~ 10~ 3 E. 



1. Gamma-Ray Emission of Pulsars 

There are now at least seven spin-powered pulsars known 
to be emitting pulsed gamma radiation. Detailed spectral 
information has been obtained by observations with the 
Compton Gamma Ray Observatory in the energy range 
of up to 10 GeV. Several of these pulsars have a break 
in the spectrum at about 1 — 2 GeV. This is also the en- 
ergy where the gamma-ray pulsars emit most of the spec- 
tral power. The gamma ray luminosity is correlated to the 
square root of the spin-down power of these pulsars (with 



some scatter). For specific informations about the proper- 
ties of individual gamma ray pulsars we refer to Thompson 
et al. (1996), Fierro et al. (1998), Thompson et al. (1999), 
and Kanbach (2001). 

The models proposed for high-energy pulsars can be 
devided into polar cap models (e.g., Daugherty and Hard- 
ing 1996; Sturner and Dermer 1994, Usov and Melrose 
1996) and outer gap models (e.g., Romani and Yadi- 
garoglu 1995; Romani 1996). In the first group the ac- 
celeration of particles and the gamma ray production is 
done in the open field line region above the magnetic pole 
of the neutron star. In the second group this takes place 
in the vacuum gaps between the neutral line and the last 
open field line in the outer magnctosphcre, near the light 
cylinder. Depending on the distribution of the electric field 
parallel to the magnetic field along and across the gap, a 
power-law energy distribution of primary particles can be 
expected that emit gamma radiation. Here we discuss the 
effects of subsequent cooling by curvature radiation. 

1.1. Primary Particles, and Particle Energy Spectrum 

Particles entering the starvation gap will be accelerated 
by the strong electric field along the magnetic field. The 
accelerating electric field given in Cheng, Ho, and Ruder- 
man (1986) can be rewritten as 
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where 8 ~ 0.3 is the ratio of the gap thickness to the light 
cylinder radius. The particles gain energy at a rate P ga in = 
ceEo. The energy of these primary particles is limited by 
curvature radiation losses, that depend on the bending of 
the field lines. The equilibrium between gains and losses 
determines the maximum energy. The curvature radius R c 
of the last open field line is of the order of the light cylinder 
radius Ri c = c/Q, with SI = 2it/P, and increases towards 
the neutral line. This becomes even more important for 
a (nearly) perpendicular rotator. By setting R c = kR\ Cl 
with k ~ 1, the power emitted by a single particle is then 
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The typical frequency emitted is 



4ttP c ' 4tt ' 



(3) 



Equating P ga i n with Pi oss = P cr , this results in an equilib- 
rium Lorentz factor for the primaries of 



7max = 4.6xl0 7 (fc*)V2«-3/4 



1/4 



(4) 



where k ~ 0.5 is the ratio of the inner starting radius of 
the gap to the light cylinder radius. 

The emission of gamma ray photons by the primaries 
leads to a pair creation cascade. These pairs, if produced 
inside the gap, will also be accelerated, start to shield the 
electric potential in the gap. As their density increases 
and becomes of the order of the Goldreich- Julian density 
tiqj = QB/2irec this finally closes the gap. The outward 
moving particles are no longer accelerated and only loose 
energy via curvature radiation. 

The pairs generated outside the gap with a finite pitch 
angle due to aberration effects will emit synchrotron ra- 
diation that can be seen in the MeV, X-ray and optical 
bands (Romani, 1996; Crusius-Watzel, Kunzl, and Lesch, 
2001). The electric field available for the acceleration of 
particles depends on the details of the cascade process, 
which leads to a quasi stationary state where the electric 
field is a function of the location along and across the 
gap. Although attempts have been made to understand 
the cascade and outer gap closure, no solid model exists. 
Therefore we just assume a power-law energy distribution 
of the accelerated particles, 



Q{i) K 7 S , 7min < 7 < 7 



(•5) 



with an index s, that must be determined from the obser- 
vations, rather than from theory in the light of the con- 
siderations given above. Since v oc 7 3 the ratio 7max/7min 
needs to be only about 5 to cover a frequency range of 
more that two decades. 

1.2. Curvature Radiation Cooling 

The timescale for curvature cooling of relativistic electrons 
or positrons near the light cylinder of the pulsar magne- 
tosphere is given by 



r = - = 4.1 X 10 
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where 7 = P cr /mc 2 and equation (||) has been used. We 
estimate the Lorentz factor 7b above which we expect a 
spectral break by setting the cooling timescale equal to 
the 1 time needed for the particle to move out of the mag- 
netosphere into the wind region, t csc = fR\ c /c = f/Q, 
where / ~ 0.1 is the fraction of the light cylinder radius, 
corresponding to the outer gap region in which the parti- 
cles are not accelerated anymore but only radiate. From 



this consideration we find that cooling will be important 
for particles with 7 > 7b, where 

The frequency at which the spectral break occurs is then 



m 3 _ 9mc 3 k 



(8) 



which corresponds to a frequency and photon energy of 
v h = 3.8 x 10 23 (fc//_i)Hz and 



E h = hvb = 1-6 GeV 



k_ 
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respectively, depending only on the geometrical factors k 
and /. In particular, the break energy is independent 
of the magnetic field strength and period. 

1.3. Gamma- Ray Spectrum 

For frequencies below the break the spectrum can be 
calculated from the particle energy distribution accord- 
ing to equation (^). By using the monochromatic ap- 
proximation for the spectral power of a single lepton, 
Pv(l) — Pcr{l)5[ v — v c{l)]i it follows that 



I v cc / d 7 Q(7)^(7) 
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where 70 oc u 1 ^ 3 is the Lorentz factor at which the argu- 
ment of the (^-function becomes zero. Wc then find that 



(11) 



This corresponds to a photon energy index of <p(u < z/b) = 
— (s+l)/3, where N-^ oc £T. The Crab, Geminga, and Vela 
pulsar have photon indices -2, -1.4, and -1.6, respectively, 
below the spectral break (Fierro, et al., 1998). This then 
gives the particle energy indices s of the primaries below 
the break energy as 5, 3.2, and 3.8, respectively. 

For frequencies above the break we first have to find 
the distribution function modified by radiation losses. 
Assuming that a spectrum (3(7) is injected at a rate 
0(7) = Q( 7 )Aesci the cooled spectrum is given by the 
following expression 



diQ(i) 



(12) 



from which wc find Af( 7 ) oc 7 ( s + 3 )^ gy usm g 
monochromatic approximation again, wc finally get 



I v OC J/~( s+1 )/3 („ > y h ) , 



(13) 
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corresponding to a photon index of (f>(y > z/t>) = — (s + 
4)/3. Therefore the spectral index of the gamma radiation 
changes by 



-1 



(14) 



at the break frequency (eqn. ^) in the direction of increas- 
ing frequencies. 

The pulsar PSR B1706-44 shows exactly this behavior 
(Thompson et al., 1996). The gamma ray photon spectrum 
changes by A<fi = —1, from <f> = —1.25 below 1 GeV to 
cf> = —2.25 above 1 GeV, with the spectrum extending for 
more than a decade in energy, both, below and above the 
break. From this we find s = 2.75 for this pulsar. 

The emitted power per logarithmic frequency band 
scales as vl v cx j/ 5 ~ s )/ 3 below the break and vl v oc 
j/ 2 - s )/3 a bove the break. For 2 < s < 5 the gamma ray 
energy distribution thus is rising below and is falling above 
the break energy. The maximum power is then emitted at 
the break energy. 

The highest frequencies are produced by the primary 
particles that feel the vacuum potential. The spectrum of 
the radiation emitted by them cuts off exponentially above 
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1.5 x lO^S^k^K-^B^PZ^TLz 



The change in the spectrum is dominated by cooling only 
if the exponential cutoff is at higher energies than the 
break, v max > v^. The cooled spectrum is nicely seen in 
PSR B1706-44, where the cutoff must be above 10 GeV. 
The situation is less clear for the other gamma ray pulsars, 
where the energy range between the break and the cutoff 
is either smaller or not existent. It can be expected that 
with the GLAST mission more pulsars with a clear break 
will be found and also that a smaller range of the cooled 
spectrum can be detected. 
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Fig. 1. The squares show the observed high energy lumi- 
nosities versus the parameter combination Bo/P 5 ^ 3 . The 
line represents the predicted luminosity at the break fre- 
quency (equation^), i.e., when radiative cooling is impor- 
tant. 

to B 2 /P 4 , the gamma luminosity thus scales according to 
L 1 oc £/i/2pV3. The observed high -energy luminosities of 
the seven gamma ray pulsars indicate a trend L 1 oc E 1 / 2 , 
and the small additional dependence on P 1 / 3 is still m 
agreement with this correlation. If we instead insert the 
maximum Lorentz factor (Q) into the power formula, the 
result is 



L 1 = 7.9 x 10 



34 
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1.4. Luminosity Estimates for the Gamma- Ray Emission 

We now estimate the power L 1 emitted by pulsars in the 
gamma-ray region at the break frequency very roughly as 
the product of the volume, the Goldreich-Julian number 
density near the light cylinder uqj oc Bq/P 4 , and the 
radiation power of a single particle P cr oc 7 / P 2 , 



L 7 = WlG J -Per , 



(16) 



where the volume is assumed to be a shell near the light 

?3 



cylinder with V = Att/R^ oc P 3 . From this we find 



L-, 
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where the break Lorentz factor from equation (R) has 
been used in the curvature radiation power (^J). Since 
the spin-down luminosity E of a pulsar is proportional 



and is thus linearly dependent on E. 

In Figures (1) and (2) we compare the model results 
with the data points. The latter have been taken from 
Thompson ct al. (1999). These observed luminosities are 
integrated fluxes above 1 eV. It is clearly seen in Figure (1) 
that if the maximum of the luminosity is determined by 
the break, a reasonable correlation is found, not only for 
the slope but also for the absolute values. However, the ob- 
served luminosity of Geminga is smaller by about a factor 
of 5, which might be due to the fact that this pulsar emits 
20% of its spin-down power in gamma rays. This means 
that the model curve would give an efficiency of 100%, 
which then causes a significant nonlinear back-reaction 
on the rotation of the pulsar and on the formation of 
outer gaps, with the effect that the gamma luminosity of 
Geminga is reduced below the expectation from the lin- 
ear model. Another point worth to note is that the Crab 
pulsar emits it maximum spectral power at ~ 100 kcV. 
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Fig. 2. Here the data points are plotted versus Bq/P 4 . 
The solid line is given by L 1 cx E 1 / 2 cx Bq/P 2 and the 
dotted line represents the predicted luminosity at the max- 
imum frequency (equation p"5| ), i.e., when radiative cooling 
is not important. 



Therefore the data point of the Crab pulsar that gives the 
luminosity in the GeV range is actually somewhat lower 
than plotted here. With these remarks in mind our simple 
model is in good agreement with the data. In Figure (2) 
the luminosity is plotted versus Bq/P 4 . The solid line rep- 
resents the dependence logL 7 = 20 + 0.5 log(£>Q/P 4 ), and 
the dotted line is the relation found, if the peak luminosity 
would be reached at the cutoff energy. The latter clearly 
overpredicts the observed values, especially for the high B 
- low P pulsars, since the slope of the curve is too steep. 

2. Synchrotron Model for the Infrared, Optical 
and X-Ray Emission of the Crab pulsar 

The Crab pulsar has a continuous spectrum from the op- 
tical to X-rays and 7-rays with different power-laws (Lyne 
& Graham-Smith 1990). The spectral index a, defined as 
I v cx v~ a , varies from a = 1.1 in the 7-rcgion via 0.7 
in the hard X-ray region to 0.5 at soft X-rays (Toor & 
Seward 1977) and zero at optical frequencies (Percival et 
al. 1993). In the far infrared region the spectrum is in- 
verted a = — 2 and cuts off sharply towards lower fre- 
quencies. This behavior is not accompanied by dramatic 
pulse profile changes, as one would expect from saturation 
or self-absorption effects. Self-absorption should first influ- 
ence the peak intensity. The infrared spectrum measured 
by Middlcditch, Pennypacker, & Burns (1983) is still not 
explained. In the following it is shown that the theory of 



optically thin synchrotron radiation at very small pitch 
angles gives a possible solution to this problem. 

2.1. Synchrotron Emission at Small Pitch Angles 

The emissivity at very low pitch angles ("J < 1/7) of a 
single particle is given by (Epstein 1973) 



e*(0,7) 
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with v B = eB/2-n-mc = 2.80 x 10 12 5 6 Hz. This formula 
has to be applied when angles 9 < 1/7 are resolved in 
the observations, i.e. when the pulse width Ac/) becomes 
comparable to the emission cone angle of the particles, 
I/7 > 2wA<p. Since the maximum of the emission is in 
the forward direction (9 = 0) the emission is dominated 
by the field lines that point towards the observer at each 
phase of the pulse. When the typical angle of emission 9 is 
small compared to the pulse width it is useful to integrate 
the emissivity over all angles (averaging): 



£1/(7) = 2-7T / e u sm9 d9 
27r 2 e 2 * 2 
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2.2. Infrared Spectrum 



The emission from a power law distribution of relativis- 
tic particles, N(-f) = 7Vo7 _s , m the direction along the 
magnetic field 9 = is given by 



I v cx / e„(0,7)iV(7)d7 



so that 
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(22) 



The spectrum rises very steeply, in fact I v cx v 2 for s = 2. 
The Lorentz factors in the case of the Crab then have to 
be of the order of 7 ~ l/27rA</> ~ 20, which is is consistent 
with the requirement that 2~/i>b lies in the infrared part 
of the spectrum. 

2.3. Optical Spectrum 

In case of somewhat higher Lorentz factors the single par- 
ticle beam angle becomes smaller than the phase angle of 
the pulse and the total (or angle averaged) power spec- 
trum has to be applied. For simplicity the monochromatic 
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approximation is used: 

dE „, 

e v ^-—5{ v -2 1VB ). (23) 

This is a good approximation since the emission is sharply 
peaked at 2^vb and one may just put all the emission at 
this frequency. The result for the optical emission is then 

I v «Je v N{ l)dl = ^-{^) • (24) 

In the case of s = 2 this gives the observed flat, a 0, 
optical/UV spectrum of the Crab pulsar, with 7 ~ I0 2 . 

2.4- Soft X-ray Spectrum 

Soft X-ray emission in this model comes from the particles 
with larger pitch angles ;§> 1 /j) , radiating at a typical 
frequency v c = |^ v f'7 2 The spectrum in this large pitch 
angle limit is given by 1„ oc ;/~(* -1 '/ 2 . This yields the 
observed l v oc v~ ' 5 at soft X-rays for s = 2. 

2. 5. Luminosity Estimates for the Optical and X-Ray 
Emission 

The energy radiated by the pulsar in the optical and in 
the X-ray band can be estimated simply as the number 
of particles times the energy radiated by one particle L = 
MriGjVP S y n , where M ~ jp/j is the multiplicity of the 
relativistic particle density within the volume V as com- 
pared to the Goldrcich- Julian density ncj- The volume 
is estimated as a spherical shell within a fraction / of the 
light cylinder radius V = 47ri? 2 c -/i?i c = 1.4x 10 30 /P 3 cm 3 . 
The frequency of the optical emission (at small pitch an- 
gles) is given by 

v op t = 5.6 x 10 14 7opt, 2 5 6 Hz (25) 

and the optical luminosity: 

iopt = lO 33 ^!^!*,; 1 _ 3 *opt,-2^opt,14.7 

x bVLbP-m erg' 8 ' 1 ■ ( 26 ) 
The frequency of the X-ray emission is 

= 2.8 x 10 17 ^ x ,_ l7 2 , 3 B 6 Hz (27) 

and the soft X-ray luminosity evaluates to 

L x = 6 x lO 35 ^/-!^^-!^ 

x Bo.ia.s^.s ergs' 1 . (28) 

This can now be compared with the spin-down luminosity 
I/ S d of the Crab pulsar due to magnetic dipole radiation: 

L x w 2 x 10- 3 aJ_ 1 /_i*- 1 _ 3 *^ 1 ^ / 1 2 7 E ■ (29) 

Since this relation is now independent of Bq and P, it 
should apply to any pulsar. Becker & Trumper (1997) have 
found such a relation observationally for X-ray selected 
pulsars. More details of the model presented here can be 
found in Crusius-Watzel, Kunzl. & Lesch (2001). 



3. Discussion and Conclusions 

We have found that curvature radiation cooling results in 
a universal spectral break energy in the photon spectrum 
at ~ l GeV in gamma ray pulsars. This break depends 
only on geometric factors, but does not depend on other 
pulsar parameters as the magnetic field strength or the 
rotation period. The maximum power is emitted at this 
energy. In fact in several of the gamma-ray pulsars a spec- 
tral break or cut-off can be observed at this energy. A 
photon spectral index change of Ac/> = — 1 is expected, 
when going from lower to higher energies. This is clearly 
seen in PSR B1706-44 (Thompson et al., 1996) and we 
expect that during the GLAST mission more pulsars with 
this spectral feature will be detected. 

The high energy luminosity calculated with cooling is 
proportional to Bq/P 5 ^ 3 , which compares well with the 
observed values. Romani (1996) has obtained cutoff en- 
ergies from a curvature radiation model without cooling 
that also lie in the observed range, but they depend on 
the pulsar parameters Bq and P, since they are connected 
to the maximum energy of the primaries. Using the cut- 
off particle energies for the luminosity estimate we find 
that the power output of pulsars in high energy photons 
scales as oc _B 2 /P 4 , i.e. as a constant fraction of the spin- 
down power E, which ovcrpredicts the 7-ray luminosity 
for young pulsars. But note that in deriving this result it 
was assumed that the density of the emitting particles is of 
the order of the Goldrcich- Julian density, which might not 
be true in general. Harding (1981) estimated the gamma 
ray luminosity from a polar cap model by assuming that 
the power put into particles is completely converted into 
radiation and found a scaling oc BqP~ 2 . In the outer gap 
model with a power-law energy distribution of relativistic 
particles (due to gap closure), the "luminosity" in parti- 
cles is dominated by the density at the lowest energy. On 
the other hand the maximum spectral power is dominated 
by the high energy particles with 7 = 7b or 7 = 7 max - So 
it is not immediately clear what fraction of the particle 
luminosity is converted into radiation. 

According to Kanbach (2001) the observed luminosity 
trend L 7 oc E 1 / 2 can be caused either by a proportionality 
to the number of particles, which would all be accelerated 
to the same energy, or by the same particle flow from dif- 
ferent pulsars with correspondingly different acceleration 
energies. In our opinion the reality is somewhere in be- 
tween, depending on the details of the pair cascade and 
the gap closure. Only when the energy (depending on the 
locally available electric field and curvature radius) and 
the density of the relativistic particles, including its spa- 
cial dependence is known, we can get more insight into 
the spectrum and the luminosity of gamma-ray pulsars. 

Synchrotron emission of particles with small pitch an- 
gles with a single power-law energy distribution, N(j) oc 
7 -2 reproduces the key properties of the Crab pulsar spec- 
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trum and the luminosity from the infrared up to soft X- 
rays. 
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